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I""! Abstract. Functional renormalization group (FRG) is applied to the three-body 

5^ scattering problem in the two-component fermionic system with an attractive contact 

biO interaction. We establish a new and correct flow equation on the basis of FRG and 

show that our flow equation is consistent with integral equations obtained from the 

Dyson-Schwinger equation. In particular, the relation of our flow equation and the 

^ Skornyakov and Ter-Martirosyan equation for the atom-dimcr scattering is made clear. 
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1. Introduction 

Few-body physics often provides fruitful information to study many-body physics in 
non-relativistic systems. A typical example is the two-component fermionic system 
with an attractive contact interaction. At sufficiently low temperatures, the superfluid 
phase of this system shows crossover from the Bardeen-Cooper-Schrieffer (BCS) type 
superfluidity to the Bose-Einstein condensation (BEC) as the attraction becomes 
stronger [lf[3]. In the BEC region, existence of composite bosons, or dimers, implies 
that low-energy excitations must be described by those dimers. Since dimers dominate 
the system in the BEC region, the scattering property between two dimers provides 
important information for quantitative understanding of BEC (see e.g. [4]-[6]). 

Functional renormalization group (FRG) (7]-[9], which has been recently developed 
as a non-perturbative method of quantum field theory, provides a unified description 
between few-body and many-body physics. Indeed, it has been applied not only to study 
few-body physics in many different contexts such as three-body or atom-dimer scattering 
problems 10,11 , trion formation [l2|, and Efimov physics 13 15 , but also to study 



many-body physics, such as BCS-BEC crossover 16 17 . FRG is a formalism to calculate 
one-particle-irreducible (IPI) effective action T[iIj] of fields ijj,ip in a non-perturbative 
way. In this formalism, we define the IPI effective action Tklip] by adding the two-point 
function ipRkip to the classical action S[ip], where k is a parameter controlling the scale. 
Then the functional Tk[ip] obeys the exact one-loop flow equation (^-[9) 

1 



dkVklij] = -STr 



^(2)r 



dkRk{T'j:'[^] + Rkr , (1) 



where F^ [ip] is the second derivative with respect to the fleld ip: 

This equation is often called the Wetterich equation and holds an exact one-loop 
property. If the regulator term Rk vanishes as A; — > 0, then F^ reduces to the usual 
IPI effective action F, which contains all the information of quantum fluctuations. 

FRG must be consistent with the Dyson-Schwinger equation, since both methods 
provide exact relations among Green functions. Theoretically, we can prove that the 
Dyson-Schwinger equation can be regarded as an integrated equation of the Wetterich 



equation 18 . However, in previous studies 10 15 on few-body physics with FRG, 
flow equations in general contradict Dyson-Schwinger equations. Those flow equations 
were constructed based on observations that the Wetterich equation has exact one-loop 
structure and particle-hole loop vanishes in the vacuum for non-relativistic physics. 
However, some one-loop diagrams in the flow equation, which look like particle-hole 
loops at flrst sight, do not vanish even for non-relativistic scattering problems in the 
vacuum. As a result, a flow equation for 2r2-point IPI vertex depends on {2n + 2)-point 
IPI vertex. 

In this study, we consider three-body problems for a contact and attractive 
interaction of two-component fermions, and establish the correct flow equation of FRG 
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in a concrete way. In order to construct the correct and closed flow equation, we 
combine diagrammatic techniques with FRG so as to identify feedbacks from higher 
point vertices. Our new flow equation is a closed equation and can be shown to 
be equivalent to the integral equation for three-body problem originally derived by 



Skornyakov and Ter-Martirosyan 19 



Outline of this paper is as follows. In section [2| we fix our notation and briefly 
review the result on two-body scattering problem with an attractive contact interaction 
of two-component fermions. In section [3} we consider its three-body scattering problem 
using FRG and analyze detailed structure of the flow equation. There we will find that 
feedback from higher point vertices do not necessarily vanish, however diagrammatic 
considerations will open a way to solve this hierarchy problem. As a result, we establish 
a new and correct flow equation for three-body problems in non-relativistic few-body 
physics. In section |4| we discuss the relation between our flow equation derived in section 
[3] and the integral equation for three-body problem originally derived by Skornyakov 



and Ter-Martirosyan 19 . This makes clear the connection between few-body physics in 
FRG and other conventional methods for scattering problems. Finally, we summarize 
our result and comment some perspectives in section [5] 

2. Two-body scattering problem 

Before going into three-body scattering problem, we need to consider and solve two- 
body scattering physics in the vacuum. We consider two-component fermionic systems 
with a contact interaction, where the bare classical action of that model is given by 



S= d^x 



2 



i^{x) [^r- — - 1^] i^{x) + gi)^{x)^i)^{x)i)i{x)i)^{x) 



(3) 



In tej), '?/' represents a two-component fermionic field, ip its conjugate, r the imaginary 
time, m the mass of particles, and yU the chemical potential. In order to treat scattering 
physics in this way, we need to put the system in the vacuum (the temperature T = 0, 
and the particle number density p = 0). 

In the effective action F [■?/'], all possible terms allowed by symmetries can appear in 
general. We introduce vertex functions as coefficients of the expansion in terms of the 
fields: 

oo 

r[^] = E M^ / (2^)'^'(Pi + ■ ■ ■ +Pn - gi qn) 

^^l l^-J Jpi,...,p„;qi,...,q„ 

X r;^r.:'Jr (Pl' • • • ,Pn; gn, • • • , gi)'^p,,a^ ■ ■ ■ V'p„,a„V'<7n,/3n " " " V'giA' (4) 



where / = / d'^p/(27r)'^, ipp^ipp are Fourier components of V^(x),'?/'(a;), and aj,/3j denote 
spin indices. Especially the four-point vertex rjj consists of a spin-singlet part and 
a spin-triplet part due to spin-5'f/(2) symmetry of the classical action pi). Since the 
contact interaction does not have relative momentum dependence, only the spin-singlet 
part rf (P) survives in the four-point effective couplings, where the superscript 5* denotes 
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-dkVliP) = dk 




Figure 1. Diagrammatic expression of tlie flow equation of fom'-point couplings. Each 
line represents fermionic propagator and the gray square vertices represent negative of 
four-point vertices — Ff (P). 
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Figure 2. Identification of the four-point coupling r^{P) and the bound-state 
propagator D{P). This diagrammatic identification will be used later for convenience 
even if the bound state does not exist. 



the symmetric part in terms of relative momenta and the argument P denotes the center- 
of-mass momentum. 

The flow equation for the four-point coupling Tf (P) is 

Here we introduced the /c-derivative dk which acts only on an explicit fc-dependence 
of Rk in propagators and Gq^{1) = il^ + P/2m — ft is the free inverse propagator. 
Diagrammatic expression of (tsl) is given in Figure [I] 

Since the four-point couplings does not depend on relative momenta at the 
bare scale, relative momentum dependence of Ff does not appear at any scale k. 
Renormalization condition of the four-point vertex is given by 

F^(P° = -2ifx, P = 0) = ATias/m, (6) 

where as is the s-wave scattering length. For an attractive contact interaction, the 
condition as < implies the absence of bound states and we can put /x = in realizing 
p = at T = 0. If a^ > 0, there exists a bound state with the binding energy l/ma| 
and we should put /x = — l/2ma| to realize p = and T = 0. 

If the bound state exits, we can identify the boson propagator by looking at 
its pole in the effective vertex function. Then, we can regard the fermion four- 
point vertex as in Figure ^ In the left hand side the vertex represents ^JStx /m?as 
and the arrowed double line represents the rest of Ff(P). If a^ > 0, a bosonic 
bound state exists, and the double line at /c = becomes the dimer propagator 

D{P) = - (l + ^l + maKzPO + H)') /[2{iP^ + P^/Am)]. 

The flow equation of Ff, ([s]), and graphical notations introduced in Figures [I] and 
[2] will be used in following sections. 
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Figure 3. Flow equation of the six-point IPI vertex. Each square vertex represent 
negative of corresponding IPI vertex function. Here we should regard that external 
legs of each diagram are antisymmetrized. 







Figure 4. Decomposition of the six-point vertex FfcTIf by extracting external two- 
body scattering part Ff . 



3. Three-body scattering problems in FRG 

Let us go into three-body scattering problems. In case the bosonic bound state exists, 
calculating six-point IPl vertices is necessary for solving atom-dimer scattering problem 
in purely fermionic framework. The following discussion works well even when dimer 
does not exist, but just for clarity we consider the case where dimer exits. 

Let us start with the flow equation of a six-point IPI vertex LfcTIi. The 
diagrammatic expression of its flow equation can be found in Figure |3} and its analytic 



expression will be given in Appendix A There are two other diagrams in the right hand 
side of Figure p^ when we expand the Wetterich equation (llj) in terms of fields ip^ ip, but 
diagrammatic consideration shows that they vanish because they necessarily contain 
particle-hole loops. The first diagram in the right hand side of Figure |3] may seem to 
vanish at first sight since it also looks like a particle-hole loop, but it can contribute due 
to non-locality of effective vertices. We will find in the next section that the fiow equation 
for the three-body problem can be consistent with the Dyson- Schwinger equation only 
when this term is taken into account. 

If we use the fiow equation in Figure |3| it depends on eight-point IPI vertices and 
the fiow equation of eight-point vertices again depends on ten-point vertices. In this 
way, the hierarchy of the Wetterich equation ([I]) appears even if we would like to solve 
non-relativistic few-body problems. On the other hand, as a general property of non- 
relativistic few-body physics, we can solve n-body problems without considering (n+l)- 
body problems. These two statements can be consistent, and the diagram containing 
eight-point IPI vertices in Figure |3] turns out to be written down within three-body 
scattering matrix elements as we will see. 

It is convenient to introduce the decomposition of the six-point IPI vertex Pfe j|+ 
by amputating external two-body T-matrices. Diagrammatic representation of this 
decomposition is given in Figure |4| and according to this decomposition we define a new 
vertex function Afc so that 

^kVniPuP2,P3;P3,P2^P'l) (7) 
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+ 




Figure 5. All possible diagrams for the three-body sector in eight-point IPI vertices. 
In these graphs, the arrows in double lines are omitted just for simplicity of diagrams, 
but there is no room for any ambiguity. 



dkRk 



Figure 6. Example for getting the feedback terms from diagrams in Figure [5J 





= -rf (P2+3)Afc(pi,P2+3;P2+3.K)rf(P2+3) + rf(pi+3)Afe(p2,pi+3;P2+3.Pi)rf (P2+3) 

+rfc(p2+3)Afc(pi,P2-f3;Pl+3,P2)rf(K+3)-rf(Pl+3)Afe(p2,Pl+3;PU3>P2)rf(Pi+3)- 

Indeed, we can show that this decomposition is possible directly from a formal 
diagrammatic expansion of six-point vertex in terms of the scattering length as- When 
dimer exists, this vertex A^ represents scattering processes between atom and dimer. 

With this decomposition, we can relatively easily represent the three-body sector 
in eight-point vertices, which can contribute the first diagram in the right hand side of 
Figure |3j All possible IPI diagrams in the eight-point vertex, which can contribute as 
a feedback, are listed in Figure |5| Notice that all these eight diagrams in Figure [5] are 
one-particle irreducible, with eight external lines, and belonging into three-body sector. 

The diagrammatic expression of eight-point vertices in Figure [5] gives an expression 
for the feedback term in the flow equation of A^. In order to obtain the feedback terms, 
we just have to do the following procedures, which are similar to those in calculating 
the amputated diagrams (see also Figure [6]): (i) Amputate the external bosonic lines in 
each graph, and (ii) close two fermion lines with the two-point function dkRk so as not 
to make particle-hole loops. This procedure is clearly possible in each graph in Figure 
[5| and then we here explicitly show that feedback from higher-order vertices is possible 
even in the vacuum. Furthermore, this procedure closes the hierarchy of the Wetterich 
equation ([I|, and we need not consider four-body problems in discussing three-body 
scattering physics. 

Diagrammatically, the flow equation for A^ can be represented as in Figure [7} "9 
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dk\f =dki V^ +£^ + rf) +9 feedback terms, 



X 




Figure 7. Flow equation of Afc. Here 9 feedback terms are not explicitly given but 
they can be found from Figure [5) 



t-^^>^' t-^^)^ 




Figure 8. Dyson-Schwinger equations for the six-point IPI vertex. Black dots 
represent the bare contact coupling. 



feedback terms" in Figure [7] are immediately obtainable from Figure [5] by following the 
procedure mentioned above. Using the momentum conservation, we can reduce the 
number of variables in A^. Due to the momentum conservation, P := pi + P2+3 = 
Pi +P2+3' ^^"-^ then for simplicity of notations we introduce a new notation 

h{q, q') ■■= h{P + q, -q; -q' , P + q')- (8) 

Then, we can write down the analytic expression of the flow equation in Figure [7| as 



dkXJq, q) = 9 feedback terms + dk [- , ^''^ !} r (- , Mfl^^ 

^ ^ Ji[G,' + Rk]{l + P)\[Go' + Rk]H-P 



q) 






[Go ^ + Rk]{-1 -P-q') [Go ' + Rk]{-1 -P- q)[G^^ + Rk]{-1 -P-q') 
In this way, we establish a closed equation for the atom-dimer scattering problem. The 
resultant equation (JQJ) no longer has one-loop structure. For the full expression of the 
flow equation, see Appendix B 



In case of the two-body scattering problem with a contact interaction, such 
feedbacks are absent. This can be understood again most easily from diagrammatic 
considerations. The diagrams in six-point vertices do not have any structures like 
those in ([s]), because such six-point diagrams cannot be IPI. Absence of the self-energy 
correction can also be understood from similar reasoning. 

In the following, we discuss properties of equation ^ and find its formal solutions. 

4. Formal solutions of three-body scattering problems 

In order to find a formal solution for the atom-dimer scattering vertex A^, let us consider 
Dyson-Schwinger equations for the six-point vertex functions. We should point out that 
the Dyson-Schwinger equation must hold even for theories with an IR regulator i?fc. 



and it can also be interpreted as an integrated flow of the Wetterich equation 18 
For the six-point vertex rfclli in our model, we have two Dyson-Schwinger equations 
for the scattering problem. Their diagrammatic expressions can be found in Figure [8j 
In deriving the second Dyson-Schwinger equation for the six-point vertex in Figure 
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X 




Figure 9. Diagrammatic expression of the integral equation for Afc 

A= \f , G= <t!^-*^, K 



X-H 



Figure 10. Diagrammatic definition of symbols in (111 



Sl we use the Dyson- Schwinger equation for the four-point vertex Ff and the first 
Dyson-Schwinger equation for the six-point vertex in Figure [8] This second expression 
of the Dyson-Schwinger equation is equivalent to the integral equation for three-body 



scattering problems derived by S. Weinberg in 20 



Combining these two Dyson-Schwinger equations in Figure [8] and applying the 
decomposition of the six-point IPI vertex in ([T]) or in Figure 111 we obtain the integral 
equation for the atom-dimer scattering vertex A^. Diagrammatically, A^ satisfies the 



integral equation in Figure [TTJ and its analytical expression is given by 

Afc(/,g') 



Afc(g,g') 



rf(-0 



[Go ^ + Rk] il + P)\ [Go' + Rk]H-P-q) 



+ 



1 



[Go ' + Rk]{-1 -P- g)[Go ' + Rk]{-1 -P-q'))' ^^°^ 
At A; = 0, this integral equation is equivalent to the integral equation for the atom- 
dimer scattering T-matrix, which was originally derived by G. Skornyakov and K. Ter- 



Martirosyan 19 . We should emphasize that the integral equation (10) holds at any 



values of the parameter k. Furthermore, the solution of this integral equation (10) must 
satisfy the flow equation ^ 

Let us explicitly show that the solution of this integral equation satisfies the flow 
equation (|9]). For this purpose, it is better to write down the integral equation in Figure 
|9] symbolically as 

X = G-K-X + G-K-G. (11) 

where we introduced a matrix notation so that matrix elements of A, G, and K are 
given by 

(A),,,' = A,.(g,g'), (12) 



(G),,, = [Go ^ + R 



1-1/ 



-/-P-g), 



and 



{K)y,i = (27r)V(/ - n 



rf( 



-0 



[Go'+RkW+py 



(13) 



(14) 
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In (11 ), "■" represents the multiplication as operators. Here we made their fc-dependence 



just for simplicity of notations. We can formally write down the solution of (11 ): 

A = (1 - G ■ K)-^ ■G-K-G = G-K-G + G-K-G-K-G + ---. (15) 

In order to derive the flow equation, let us take a derivative with respect to k of the 



both sides of the (15). Leibniz rule and formal resummation give 



dk\ = dkiG-K-G) 
+dk{G ■K)-{G-K-G + ---) + {G-K-G + ■■■)■ dk{K ■ G) 
+{G + G-K-G + ---)-K-dkG-K-{G + G-K-G + ---) 
+{G ■K-G + G-K-G-K-G + ---)- dtK ■{G-K-G + G-K-G-K-G + 



(16) 



Combining (15) and (16), we obtain 



dkX = dk{G ■K-G) + dk{G ■K)-X + X- dk{K ■ G) 

+ {G + X)-K-dkG-K-{G + X) + X-dkK-X. 



(17) 



Now it is clear that the equations in Figures^ and u^ and the equation (17) are equivalent 



as it should be since Dyson-Schwinger equation must hold at any values of k. 



5. Conclusion 

In this paper, we showed that some one-loop diagrams in the flow equation, which look 
like particle- hole loops at first sight, do not vanish even for non-relativistic scattering 
problems in the vacuum. If we drop those terms naively, equivalence between FRG 
and Dyson-Schwinger equation is lost. In order to construct the flow equation for those 
problems, we need to write down possible diagrams for higher-point vertices so as to 
close the flow equation. 

In order to construct a new flow equation of FRG for three-body scattering problem, 
we identified all Feynman diagrams for eight-point vertices which contribute to the flow 
equation for three-body scattering problems. Combining the flow equation of FRG 
with diagrammatic considerations, we derived a new and closed flow equation for the 
three-body scattering problem. Furthermore, we proved consistency between the flow 
equation for the three-body problems and the corresponding Dyson-Schwinger equation, 
and the relation of our method with the Skornyakov and Ter-Martirosyan equation is 
unveiled. 

Let us compare our studies with previous works on few-body physics in FRG. In 



the paper 11 , authors discussed the atom-dimer scattering with FRG by introducing 



bosonic fields via the Hubbard-Stratonovich transformation in the particle-particle 
channel. They showed that their flow equation for the atom-dimer scattering is 
consistent with the Skornyakov and Ter-Martirosyan equation only in a special limit, 
in which fermions are integrated out at first so that only dynamical degrees of freedom 
are bosonic dimers. However, there exit missing diagrams if more general situations 
are considered. As long as a general regulator for fermions is taken, fermions remain 



Flow equation of FRG for three-body scattering problems 10 

dynamical degrees of freedom and dissociation processes drawn in Figure [5] must be 
taken into account to establish the consistency with the Dyson- Schwinger equation. 
FRG has also been applied to few-body physics in order to reveal its universality, 



especially in the context of Efimov physics with the help of dimer fields 13 15 . In 



the context of renormalization group methods, Efimov physics can be interpreted as 



a limit cycle solution of the renormalization group transformations 21 . Limit cycles 
can appear only when the flow equation of A^ depends on itself quadratically, and such 
quadratic terms of A^ appears through feedbacks from higher point vertices if dimer 
fields are not introduced. Therefore, in order to discuss Efimov physics without dimer 
fields, use of our formalism would be inevitable. 

We expect that a formulation for four-body scattering problems within FRG is also 
given in a similar way, but this task would be a future work. Even though we can solve 
such problems by deriving integral equations directly from Dyson-Schwinger equation, 
formulating that problem in the context of FRG is important for revealing many-body 
properties of physical systems, such as BCS-BEC crossover. 



After completion of this work 22 , we found a related work, in which a correct 
decoupling theorem for few-body physics is presented [23]. In our work, the main 
motivation is to construct a new flow equation of FRG for three-body physics, which is 
consistent with the Skornyakov and Ter-Martirosyan equation, and it is different from 



that of 23 



Appendix A. Supplement for the flow equation of six-point vertex functions 

Analytic expression of the flow equation in Figure |3] is given by 
- dkTklW (pi , P2 , P3 ; P's , P2 ' P'l ) 

Tkl]\\{pi,P2,l,P3;p'3,l,P2,P'l) ^k\l\\{pi,P2,P3j;hp'3,P2^P'l) \ .. . 

[Go' + RkW) [Go' + RkW) ) ^ ■ ^ 

rf(p2+3)r4rr(pi>^>P2+3 - 1;p3,P2,p'i) ^Upi+3)'^kl\\{l,P2,Pi+3 - 1;p'3,P2,p'i) 
[Gq-^ + Rk]{l)[G^' + Rk]{P2+3 - I) [Go' + Rk]{l)[Go' + Rk]{pi+3 - 




^kl\}{pi,P2,P3;P2+3 - lJ,P'lM{P2+3) rfc||j(/,P2,P3;Pi+3 " ^ ^ P2 , ^M iP'l+3) 



+dk 



[Go' + Rk]{l)[Go' + Rk]{p'2+s - I) [Go' + Rkm[Go' + Rk]{pi+3 - 

1 f rf(p2+3)rf(p'2+3)rf(Pi + /) 



/ [Go' + RkW) V [Go' + Rk]{l + Pi-i')[Go' + Rk]{P2+3 - I) 

rf(pi+3)rf(p-2+3)rf(P2 + 

[Go' + Rk]{l + P2-v)[Go' + Rk]iPi+3 - 

rf(p2+3)rf(K+3)rf(Pi + 

[Go' + Rk]il+Pl-2')[Go' + Rk]iP2+3 - I) 

rf (Pi+3)rf (pU3)rf (P2 + /) 

[Go' + Rk]il+P2-2')[Go' + Rk]iPl+3 - I) 
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11 



where we have introduced the abbreviation P2+3 = P2+P3, Vi-v = Pi ^v'l^ stc. In section 
|3} we found physical meanings of each term when we consider atom-dimer scattering 
problems with the model given by (pi). The first two term in the right hand side of 



(A.l) represents dissociation of dimers, and they become two-loop diagrams when we 



represent then using atom-dimer scattering vertex Afc. Other terms in (A.l) represent 
rescattering between atom and dimer. 

Appendix B. Supplement for the flow equation of the vertex A^ 

We here give the full expression of the flow equation for the vertex function A^ based 
on the procedures given in Section [3j For simplicity of notations, we use the matrix 



notation introduced in equations (12), (13), and (14) 



Using the matrix notation, equation (|9| is written as 

dk\k = dk{G ■K-X + X-K-G + G-K-G) + 9 feedback terms. 

From the diagrams in Figure |5| 9 feedback terms are given as follows: 

9 feedback terms = G ■ {dk - dk)K ■ G + G ■ K ■ dkG ■ K ■ G 
+ G-{dk- dk)K -X + G-K-dkG-K-X 
+ X-{dk~ dk)K -G + X-K-dkG-K -G 
+ X-{dk- dk)K -X + X-K-dkG-K -X 
+ X-dkK-X. 



(B.i; 



(B.2) 



The diagrams in Figure 5 and the terms in the right hand side of equation (B.2) are set 



in the same order. Here we used the flow equation for the four-point coupling Ff , and 
then {dk — dk)K represent the self-energy correction of the dimer propagator in equation 

(B.3) 



(B.2). Substituting equation (B.2) into equation (B.I), we find that 
dkXk 



dk{G ■K-G) + dk{G ■K)-X + X- dk{K ■ G) 
+ G-K -OkG-K-G + G-K-dkG-K-X 
+ X-K -dkG-K-G + X-K-dkG-K -X 
+ X ■ dkK ■ A. 

This is the full analytical expression of equation |9] 
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